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Abstract 

\jj3 \ The complete group classification of a generalization of the Black-Scholes- 

. I Merton model is carried out by making use of the underlying equivalence and 
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2^ \ 1. Introduction 

In this paper we study the following semilinear partial differential equa- 
tion 

Ut + -d^x^u^^ + rxun, + f(u) = (1) 

which can be considered as a weakly nonlinear generalization of the celebrated 
Black-Scholes-Merton equation 

„,^.V.... + ™.-™^0 (2) 



* Corresponding author 
Email addresses: bozhkov@ime.unicamp.br (Y. Bozhkov), spawnOmath.upatras.gr 
(S. Dimas) 



that plays a remarkable role in financial Mathematics. 
Usually the terminal condition 

M(a;,T) = l, (3) 

describing the evolution of standard or "vanilla" products; the price of a 
zero-coupon bond (or of a financial option), u{x,t), which is exercised when 
t = T is considered [1, 2, 3, 4]. 

Furthermore, Eq. (2) may satisfy also other kinds of options, like the 
barrier option. A barrier option can be considered an exotic option and as 
such has features that makes it more complex than the "vanilla" option, 
[5, 6, 7]. The underlying idea is that now a barrier H{t) exists and when 
the asset price x crosses it, the barrier option u{x,t) becomes extinguished 
or comes into existence. Those two types are also known as down-and-out 
and down-and-in respectively. Often a rebate, -R(t), is paid if the option is 
extinguished. In what follows we shall consider the down-and-out type. 

In the context of the Black-Scholes-Merton equation, the barrier option 
is expressed by the conditions 

u{H{t),t) = R{t), (4a) 

u{x,T) = m.ax{x — K,0), (4b) 

where the barrier option u{x,t) satisfies Eq. (2) for x > H{t), t < T, T 
again is the terminal time where the barrier option is exercised and K is the 
strike price. A common assumption for the barrier function H is to have the 
exponential form 

H{t) = hKe-''\ (5) 

where a > and < 6 > 1 [6, p. 187]. 

Although Eq. (1) is just a mathematical abstraction, it can also be derived 
from the same stochastic argument as the Black-Scholes-Merton equation 
(2) is derived. In particular, by assuming that the wealth of the portfolio 
involves a nonlinear function of the value of an option u{x,t) instead of a 
linear function and utilizing the same boundary conditions (3), (4). Finally, 
the constants r and a"^ represent the risk-free interest rate and the variance 
of the rate of the return on u respectively. 

The main purpose of this paper is to carry out a complete group classi- 
fication of a generalized Black-Scholes-Merton equation of type (1). Recall 
that to perform a complete group classification of a differential equation (or 



a system of differential equations) involving arbitrary functions or/and pa- 
rameters, means to find the Lie point symmetry group Q for the most general 
case, and then to find specific forms of the differential equation for which Q 
can be enlarged, [8, p. 178]. Quite often, there is good physical or geometri- 
cal motivation to study such cases. Moreover for the Black-Scholes-Merton 
equation (2) the motivation is a financial one: for more than three decades 
this equation has shown its value and gain the trust of the market. But as 
banks and hedge funds relied more and more on this and its siblings equations 
they became more and more vulnerable to mistakes or over-simplifications 
in the mathematics involved in deriving them; it is a linear equation and 
as such can't follow the dynamic nature and the intricacies of the modern 
market. 

In [9, 10, 11] the Black-Scholes-Merton equation (2) as well as other 
linear evolutionary equations which appear in financial Mathematics have 
been investigated from the point of view of the S. Lie symmetry theory. Also, 
the barrier option was studied recently for Eq. (2) under the same prism [7]. 
The present work can be considered as an extension of that research. As 
such, the linear case f{u) = au + f3, to avoid a repetition of already known 
and established facts, will not be treated explicitly in the present paper. It 
will only briefly accounted for the purpose of completeness. 

Afterwards, for each nonlinear case of the obtained group classification 
we look for group invariant solutions taking also into account the boundary 
conditions (3), (4). 

Before we proceed with the group classification, we would like to remind 
some facts concerning the group analysis of differential equations. 

The symmetry analysis of differential equations is a method first devel- 
oped in the 19th century by Sophus Lie. One of the main benefits of this 
method is that by following a completely algorithmic procedure one is able to 
determine the symmetries of a differential equation or systems of differential 
equations. Grosso modo, the symmetries of a differential equation transform 
solutions of the equation to other solutions. The Lie point symmetries com- 
prise a structural property of the equation, in essence it is the DNA of an 
equation. The knowledge of the symmetries of an equation enables one to 
utilize them for a variety of purposes, from obtaining analytical solutions and 
reducing its order to finding of integrating factors and conservation laws. In 
fact, many, if not all, of the different empirical methods for solving ordinary 
differential equations (ODEs) we have learned from standard courses at the 
undergraduate level emerge from a symmetry. For instance, having at our 



disposal a Lie point symmetry of a first order ODE, we can immediately get 
explicitly an integrating factor. Furthermore, even the knowledge of a trivial 
solution of the equation can be used for creating nontrivial solutions by us- 
ing the equation's symmetries. And all these are due to the rich underlining 
algebraic structure of the Lie groups and algebras with which we give flesh 
to the symmetries of a differential equation. 

Another important characteristic of the symmetry method is that in some 
situations the symmetries of an equation may indicate that it can be trans- 
formed to a linear equation. In addition, its symmetries provide the means 
to construct the needed transformation. A strong indication for that is the 
existence of an inflnite dimensional Lie algebra, [12]. Recently, it has been 
shown, by using only the algebraic properties of the symmetries of the equa- 
tions involved, that the majority of the well-known differential equations used 
in economics are linked via an invertible transformation to the heat equation, 
[11]. Although for some of the studied equations a transformation between 
them and the heat equation was already known by other means, the authors 
of [11] used the algebraic properties of the symmetries of the equations and 
followed a straightforward and algorithmic approach. 

In fact, over the last forty years there was a considerable development in 
the mathematical analysis of partial differential equations which arise in fl- 
nancial Mathematics. However when one reads various papers devoted to the 
resolution of evolution partial differential equations which arise as, more or 
less, the flnal stage of the mathematical modeling of some flnancial process, 
one can see some ad hoc and naive procedures to bring the considered equa- 
tion under control. A viable alternative way is the employment of symmetry 
analysis of differential equations. 

Moreover, one valuable tool when considering classes of equations is the 
use of the equivalence or admissible transformations, [13, 14, 15]. Equivalence 
transformations of a class of differential equations are point transformations 
that keep this class invariant, in other words they map an equation from this 
class to another member of the same class. In the recent years equivalence 
transformations have found much application either as a stand alone analytic 
tool for the group classiflcation of differential equations, [16], or at the core 
of the enhanced group analysis, [17, 18, 19, 20]. 

One of the advantages of this approach, as already emphasized, is that it 
provides a well deflned algorithmic procedure which essentially enables one 
to find the involved linearizing transformations, conservation laws, invariant 
solutions, etc. On the other hand, the calculations involved are usually very 



difficult and extensive even for the simplest equations. Thus, it may become 
very tedious and error prone. For this reason the real progress in this area 
occurred in the last few decades with the advances in computer technology 
and the development of computer algebra systems like Mathematica, Maple, 
Reduce, etc. Based on these systems, a handful of symbolic packages for 
determining the symmetries of differential equations exists, [21, 22, 23, 24]. 
One such symbolic package, based on Mathematica, [25] , has been devised and 
developed by S. Dimas as part of his PhD thesis, [26]. The package, named 
SYM, [26, 27, 28], was developed from the ground up using the symbolic 
manipulation power of Mathematica and the artificial intelligence capabilities 
which it offers. It is being used by many researchers around the world. It 
was extensively used for all the results in the present paper, both for the 
interactive manipulation of the found symmetries and for the classification 
of the equations employing the symbolic tools provided by it. 

This paper is organized as follows. In section 2 we present the basic con- 
cepts of the Lie point symmetry approach to differential equations used in 
the paper. In section 3 we obtain the complete group classification of gener- 
alized Black-Scholes-Merton equations of type (1). In section 4 we provide 
invariant solutions for each non-linear case found by the group classification 
under the two specific boundary problem studied, the "vanilla" option and 
the barrier option. Finally, in section 5 we discuss the obtained results as 
well as possible applications. 

2. Preliminaries 

In this section we expose some notions of the modern group analysis that 
will be encountered in the main sections of the article suitably adapted to the 
article's needs. For a full treatise in the subject there is a wealth of classical 
texts that encompass all aspects of the theory, [8, 12, 29, 13, 30, 31]. 

A Lie point symmetry of Eq. (1) is a one-parameter transformation of the 
independent and dependent variables 

X = x{x,t,u, e), 

i = i{x,t,u, e), (6) 

u = u{x,t,u, e), 

that keeps (1) invariant: 

Ui + -a'^x^usx + rxuai + f{u) = 0. (7) 



By substituting in (7) the point transformations (6) and using (1) we obtain 
an equation called the symmetry condition. From it the exact form of the 
point transformations is obtained. However, trying to determine the trans- 
formation (6) from (7) is more challenging than trying to solve the Eq. (1) 
itself! 

The novelty behind S. Lie idea resides on linearizing the problem of de- 
termining the symmetries of an equation. To do that, the notion of the 
infinitesimal generator is introduced. Namely, the infinitesimal generator is 
a differential operator 



X 



^id 



d 



d 



dx dt du 



with the functions ^^ = C,^{x,t,u) and rj = ri{x,t,u), called infinitesimals, 
defined as 



^ de 



e=0 



,2^dt_ 
^ de 



e=0 



du 



e=0 



An infinitesimal generator of this type determines a Lie point symmetry of 
Eq. (1), if and only if, its action on the equation will be, modulo the equation 
itself, identically zero, that is: 



X(2) 



Ut + -a'^x'^u^^ + rxu^j, + f{u) 



0, 



(9) 



(1) 



where X*^^) is the second order prolongation of the operator X given by 






t, ij 



1,2 



(10) 



and 



.(1) 



.(2) 



.(1) 



VI ' = DiT]- {Di^^)uj, and r]^'^ = Di^r]l> - {Di^^^)ui,j, i,ik,j 



1,2 fill 



where u, = ^, Ui^i^ = q^^, i,ij = 1,2 and {x^,x^) = (a;,t), where the 
Einstein summation convention implied over the indexes. From the Eq. (9), 
called linearized symmetry condition, an overdetermined system of linear par- 
tial differential equations emerges. By solving this system, called the deter- 
mining equations, we find the infinitesimals and hence we obtain the point 
symmetries of the equation. The group classification occurs in that phase. 
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The determining equations contain also the parameters a, p and the function 
f{u). The group classification is performed by investigating each case where 
specific relations among the unknown elements remove equations from the 
set of determining equations and hence, expanding the solution space. 

This set of symmetries, represented as differential operators, form a Lie 
algebra. The system of ODEs (8) with the addition of the conditions x\e={) = 
X, i\£=o = t, u\s=o = u forms a well posed initial value problem. By solving 
it, we can obtain the corresponding local continuous transformations, which 
form a Lie group. This process is called exponentiation. Henceforth, we shall 
identify a Lie point symmetry with its infinitesimal generator. 

Having the symmetries, as a Lie algebra, there is a wealth of things that 
can be done with. In the present paper, we use them to obtain invariant or 
similarity solutions of the Eq. (1). By invariant solutions we mean solutions 
of (1) that are invariant under one of the found symmetries X, e.g. 

X[u-v{x,t)]\u=^i.,t)^^- (12) 

The Eq. (12) is a linear PDE called invariant surface condition and by solving 
it we obtain a way to reduce Eq. (1). For example, the symmetry ^ yields the 
invariant surface condition Uf = 0. Solving it, we get the invariant solution 
u{x,t) = (/){x) which, in turn, can be used to reduce the order of Eq. (1) 
effectively turning it from a PDE to an ODE. Similarly, when we look for 
a similarity solution of Eq. (1) along with a initial/boundary condition we 
have to choose the subalgebra leaving also invariant that condition and its 
boundary: 

Xit-T)\,^^^0 (13) 

and 

X{u-1)\,^^^0. (14) 

for the boundary condition (3). And 

X(x-iJ(t))U^(,)^0 (15) 

and 

Xiu-Rit))l=H(t)^^- (16) 

for the boundary condition (4a). 

For obtaining the equivalence transformations there two possible roads. 
The first one is by a direct search for the equivalence transformations, an 
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approach that gives in theory the most general equivalence set of transfor- 
mation, the equivalence group. But it has the same pitfalls as trying to obtain 
the symmetry group for the equation as briefly discussed already. And the 
second approach is through the calculation of the equivalence algebra from 
which the continuous equivalence group can be obtained. In the present work 
the second road will be realized by complementing the usual prolongation of 
the infinitesimal generator with a secondary prolongation, [29]. 

To calculate the equivalence algebra, an extension of Eq. (1) must be 
considered with the arbitrary elements a,r,f, now functions of x,t,u, and 
by including the following constraints on them, 

(Tx = (^t = (^u = r^ = rt = Tu = fx = ft = 0. 
For this extended system the infinitesimal generator is 

A- = {'^+5.9 a+^.a+^.9+^3 9 (17) 

OX ot ou oa or of 

where now the coefficients of this operator depend on the extended space: 
C = Cix,t,u,a,rJ), T] = 7]{x,t,u,a,rJ) and 0* = (j)'{x,t,u,a,rj). The 
second prolongation needed to obtain the determining equation now becomes 

'' dui '''''dui.i, ^' daj ^' drj 

+ 0?'^'^T^, h^k = l,2, J = 1,2,3 (18) 



' df: 



J 



where the coefficients r]^ \^ij^ are calculated as usual by the formula (11) 
while for the coefficients (pi with the secondary prolongation. 



,.,<., 



Dif - (Die)pi, - {D,e)Pl - (D,<,)p{o , i, J = 1, 2, 3 



where {p^,p'^,p^) = {cr,r,f), {x^,x^,x^) = {x,t,u) and 



d_ , _d_ 

dxi ^ ^^' dp' ' 



A- = 7^ + Px. 



again the Einstein summation convention is used for the index i. After that 
point we follow Lie's algorithm as usual. 



Having the equivalence algebra by exponentiation one can obtain the 
continuous part of the equivalence group. By using the method proposed 
in [30, pp. 187 c.f.], [32] one can also obtain the discrete part and hence 
retrieve the whole set of equivalence transformations permissible by this class 
of equations. 

Another useful notion is that of the additional equivalence transforma- 
tion. An additional equivalence transformation is a point transformation 
that connects inequivalent classes of differential equations. The knowledge 
of such transformations greatly helps the classification. 

3. Group classification 

In this section we proceed with the group classification of the generalized 
Black-Scholes-Merton equation (1). 

First, the best representative for the class of equations (1) is obtained 
utilizing its equivalence algebra. To do that, the continuous part of the 
equivalence group is constructed and with its help as many as possible arbi- 
trary elements are zeroed. 

Theorem 3.1. The equivalence algebra C^ of class (1) is generated by the 
following vector fields 

dt, du, xda:, dr + txd^, fdf + udu, 

X {2rt -ta^ + 2 log|a;|) 9^ + Udt - Afdf, 

x{{2r + a^)t-2 log|a;|) d^ - 2ad^. 

Proof. By applying the second order prolongation (18) of the infinitesimal 
generator (17) to the extended system 

Ut + -(^{x, t, u^x^Uxx + r{x., t, u)xux + /(x, t, u) = 0, 

CTx = CTt = CTu = Tx = rt = Vu = fx = ft = 0, 

modulo the extended system itself, we get the system of determining equa- 



tions: 



Vsf = 0, V4f = 0, e? = 0, e? = 0, e? = 0, e?/ = 0, V3u = 0, V4u = 0, 

C« = 0, ^If = 0, ^l^ = 0, r^ii = 0, r]2t = 0, r^s^ = 0, r]4t = 0, r]i.^ = 0, 

^2,; = 0, r]s^ = 0, 7^4^ = 0, ^^ = 0, r]if + f^) = 0, r]if + f^) = 0, 

e} - rx^J = 0, e}/ - rx^J; = 0, euf - rxeuf = 0, ^i - rx^l = 0, 

2e? + r^i// + /e?/ = 0, 2^} - X (2 (r + a^) ej + a^a^e^,) = 0, 

2^ -x{2{r + a') C + xa^J = 0, 

rej + eu + Viuf + fCf - Cf + rxCf = 0, 
/e) + X {-Mj + xa' {C + Vixf + fCf)) = 0, 
4x7/3 + a {Ae + X (-2/e^ + 2^^ - 4^^ + QrxC + ^xa^C + x'^a^Cx)) = 
2r]2 - 2/77i„ - 2 feu + 2r7i, + 2/^ + 2rxvi^ + 2/rxe^ 

I _2_2/„ I r,-2 \ n 



2x7/4 + 2re + 2feu - 2frxC - 2^ + 2rxe - 2rxC + 2r^x^C + 2rx^a^C 

+2x^^77, „ + 2fxVVi, - xVVL + rxVVL = 0. 



Solving the above system the equivalence algebra C^ is obtained. D 

Lemma 3.2. The continuous part of the equivalence group, Sc, consists of 
the transformations 

i=5i + 5Qt, u = 52 + 53U, f- = r + 55, a = 67a, f = —f, 



"6 



where 6i are arbitrary constants and 63,64^,6^,67 7^ 0. 



Remark 3.1. Due to the fact that the transformation for x depends also 
on the arbitrary elements of Eq. (1), Sc is also called the continuous part 
of the generalized equivalence group. If one has chosen to assume that the 
equivalence transformations for x, t, u do not depend also on the arbitrary 
elements r, a, f , i.e. ^* = ^*(x, t,u), rj = rj{x, t, u) in (17), then the continuous 
part of the usual equivalence group Sc would be obtained. 
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Readily, using tc one can find an equivalence transformation that a — ?> -\/2 
and f — )■ 0: 



CT^-2r 



^^^y"^: V*|x|^, t = t, M = M, f = 0, a = V2, / = /. (19) 
Using (19), Eq. (1) turns into the equation 

Ul^3?Ui^\]{u)=^. (20) 

Next, using the additional equivalence transformation 

a; = log|x|, t = t, u=\x\~^''^u^ (21) 

the heat equation with a nonlinear source is obtained (for clarity henceforth 
the hats are dropped) 

ut + M.. + e-^/V(e"/'M) -\u = ^- (22) 

Therefore, without any loss of generality, Eq. (22) will be classified instead. 
In addition, the terminal condition (3) is transformed to the condition 

m(x,T) =exp(-x/2) (23) 

and for the barrier option condition (4a) we have: 

Finally, the equivalence group for Eq. (22) is given by the following theorem 

Theorem 3.3. The equivalence group, S, of Eq. (22) consists of the trans- 
formations 

X = (55/3x + (/3 - (55)(55t + 2^4, i=S5t + 6i, 
where 6i are arbitrary constants, 63,65 7^ and a,/3 = ±1. 
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Proof. The process is analogous to the one for Eq. (1) with the only difference 
that now we have only one arbitrary element, the function /. In addition, 
using the process described in [32] we find the four discrete equivalence trans- 
formations 

{x, t, u, f) -^ {I3x + (/3 - l)t, t, ae-5(/3-i)(^+*)^i, «/), 

where a, /3 = ±1. Together the two sets of transformations comprise the 
usual equivalence group of transformations £. D 

Equation (22) belongs to the class 

Ut = u^^ + f{t,x,u,u^) 

that describes nonlinear heat conductivity processes. Due to the fact that 
the above class was completely classified in [33]. Apart from mentioning the 
classification equation 

I {uf (e^/\) - e-/'f (e^/\)) (^3(t) + ^) + um) 

+ (e-^/'f {e^/\) - ^) (j-^(t) - J-4(t) - ^x (4J-^(t) + xj-^'(t))) 
+ (/' (e^/'w) - I) (j'iix,t) + i«(8J-4(t) +x(4J-^(t) +xJ-2"(t))) 

+ ^u (2 j;'(t) + X (4J^'(t) + xJ^^'it))) + J^uix, t) + J-i,,(a;, t) = 

o 

no further details of the calculations involved will be showed. We proceed 
with presenting the resulting classification. 

1. For an arbitrary / the Lie point symmetries of (22) are determined by 
the infinitesimal generators: 

2. For /(C) = ~^{a + (iC){5 + \og\a + /3C|) , /3, 7 ^ 0, in addition to the 
symmetries (25) we get also the symmetries 

^3 = ^ e^---, (26a) 

^^ ^ l^,tl_ ^ hit + x)e/^u + a{l + 7(t + x)) ^_.^^, d 
7 dx (5^ du' 

12 



3. For f{() = a{( — f3Y, a 7^ 0, in addition to the symmetries (25) we 
have the symmetry 

X3 = 2ix-t)^+ 4t^ + iit-x- A)u + 4/3e-^/2)|^. (27) 

4. For completeness, we present also the additional symmetries for the 
linear cases: 

• /(C)=/3C + «, /^T^O: 

/ '^ -t/2\ 9 

d f a{t + x)e'^^'^\ d 



X4 = 2t—- + { UX + 



dx \ (3 J du 

X, = 2x- + At- + (^^-(4/3-l)(^ + /^g^N0 ^-./2g^ 
dx dt (5 du' 

Xe = 4xt— + Ae- + -{a {2{x - l)t + x^ + (1 - A(5)e) e'^'^ 
ox at p ^ ^ ' 

+/3 [f + (t - 4/3t - 2)) tu) ^, 
d 

Xoo = J'(x,t) — , 

where the smooth function T = J^{x, t) satisfies the linear equa- 
tion 

J^t + J^xx + {a - -)J^ = 
and, 
/(C) = «: 

X3 = (w + axe"^/^)--, 
ou 

X4 = At-^ + (2e^/^xu + a(t^- x(x + 2))) e-"/^-^, 
ox ^ ^ " ou 

X5 = Ax^ + 8t^ + (2e"/2tM + a{e-ix- Q)x)) e""/^^, 
ox Ot " ou 

Xe = 12xt— + 12t^— + (3 ((t - 2)t + x^) u 
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+a {t^ - Ut"^ - 3tx^ - 2x{x{x + 3) + 6)) e^^/^) ^, 

where the smooth function J-" = J^{x, t) satisfies the hnear equation 

J^t + J^x. + {a- -)-F = 0. 

Remark 3.2. Looking at the corresponding Lie algebars for the two hnear 
cases, it is evident that they are hnked to the heat equation via a point 
transformation, [11], — an additional equivalence transformation — a fact 
well established in the literature, [9, 10]. 

4. Invariant solutions 

Having obtained the complete group classification for Eq. (22), and con- 
sequently for Eq. (1), we can look for invariant solutions under the terminal 
condition (23) and the barrier option (24): For each one of the two nonlin- 
ear cases the appropriate sub algebra of symmetries also admitted by each 
problem is found using the two required conditions (13), (14) and (15), (16) 
adapted now to Eq. (22). Finally, by using the subalgebra obtained for ev- 
ery subcase that surfaced from the two conditions a similarity solution is 
constructed. 

4.I. The terminal condition 

Case 4.1.1. /(() = a(C - /3)^ a ^ 

Let the arbitrary element of the Lie algebra spanned by (25) and (27) be 
X = ciXi + C2X2 + C3X3. Using (13), (14) and (23) we obtain the conditions: 

ci = -4Tc3 

and 

(/3 - 1)C3 = 0. 

From the above conditions two specific subcases occur: 
• Subcase 4.1.1.1. ci = C3 = 0, C2 7^ 
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For this case the only symmetry that keeps invariant both Eq. (22) and 
(23) is the X2 = 28^ — udu- 

• Subcase 4.1.1.2. /3 = 1, ci = -4Tc3 

For this case the only symmetries that keeps invariant both Eq. (22) and 
(23) are 

Zi = 2d^ - udu, 

Z2 = 2{x - t)d^ + 4(t - T)dt + {{t-x- 4)m + 46-^/2)^^ 
Case 4.1.2. /(C) = -^(a + K) {5 + log|a + /3C|) , /?, 7 7^ 

Let the arbitrary element of the Lie algebra spanned by (25) and (26) be 
X = ciXi + C2X2 + C3X3 + C4X4. Using (13), (14) and (23) the conditions are 

ci = 

and 

{a + /3) (c4 + 7 (c3 + (t + x)ci)) = 0. 

From the above conditions two specific subcases occur: 

• Subcase 4.1.2.1. ci = C3 = C4 = 0, C2 7^ 

For this case the only symmetry that keeps invariant both Eq. (22) and 
(23) is the X2 = 28^ — udu- 

• Subcase 4.1.2.2. /3 = -a ^ 0, 7^0, ci = 

For this case the only symmetries that keeps invariant both Eq. (22) and 
(23) are 

Zi = 2d^ - udu, 
Z2 = (1 - e"/2„)e-t+^*9„ 

Z3 = '-e^% + 7(t + x)e-/^.-(lH-7(t + x))^_,,,,^^ 
7 7 

The common denominator for all the above cases is the symmetry 2dx — 
udu- This symmetry gives the invariant solution 

u{x,t) = C{t)e-'''^- 
Going back to the Eq. (1) we see that it corresponds to the trivial assumption 

u{x,t) = C{t), 
i.e. solutions that do not have dependency on the variable x. 

15 



4-2. The condition for barrier option 
Case 4.2.1. /(() = a(C - /3)^ a ^ 

Let the arbitrary element of the Lie algebra spanned by (25) and (27), 
X = ciXi + C2X2 + C3X3. Using (15) and the equivalence transformations the 
first condition turns to the ODE: 

(2V2r (ci + 2tc3) + a((2- V2(t) ci + 4c2 - 2V2tac3t)) H 

+ AV2c:iH\og\H\ - 2V2 {ci + Atcs) H' = (28) 

two cases are discerned, C3 7^ and C3 = 0: 
• Subcase 4.2.1.1. C3 7^ 
The solution of (28) is 

H{t) = e-^^+(-i-')(*+-4^^) (29) 

where k, = j^, A = '^^^^'^^ and A the constant of integration. Using this 
solution, conditions (16) and (24) we get the ODE 

f3-R-{K + t)R' = 0. (30) 

Its solution is 

«W^^, (31) 

where B the constant of integration. 

Having found the functions H, R admitted by the symmetries we proceed 
with the reduction of the Eq. (22). From the invariant surface condition (12) 
the invariant solution is 



e-l(-+^)Ue^pt + F(^^)) 



uM = 4(. + t) (^2) 

Substituting (32) to (22) for this particular case for / we arrive at the re- 
duction 

lQe^^[5K{[3K + -)- 8(^ + /3«:)F(C) + e-^^F(C)' - -(F' + -F" = 0, 
a la a a 
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x+t+\ 

found in a closed form, it has the special solution 



where ( = ^^xf • Although the general solution of this equation cannot be 



F{C) = 4e^^ { /3k - ^ 



2«Cl 
Using it in combination with (32) we arrive to the invariant solution 

{t + X + A)^ 

Finally, using the equivalence transformations and the boundary condition 
(4a) we arrive to the similarity solution for Eq. (1) with f{() = a{( — f3)'^ 

u[x,tj = p — 



a (2A(7 + 72(^2 _ 2r)t + 2v^log|a;|) 

with 

H{t) = e-^^+(»-i-')(*+-4^^) 

and 

12^2 

R{t) = f3 2 

^^ A^a{t + k) {a^ - 2rf 

• Subcase 4.2.1.2. Cg = 

For C3 = the solution of (28) becomes 

H{t) = ^e^(2r-a=^)t+At (33) 

where A = r- "^^V -^ and A ^ the constant of integration. Using this 
solution, conditions (16) and (24) we get the ODE 

CiR' = 0. (34) 

Hence 

R{t) = B, (35) 

where B is constant. 



^ci 7^ otherwise C2 must also be zero. 
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Having found the functions iJ, R admitted by the symmetries we proceed 
with the reduction of the Eq. (22). From the invariant surface condition (12) 
the invariant solution is 



--^--K"^^^^) 



u{x, t) = e-'^'^F ^ / '- . (36) 



Substituting (36) to (22) for this particular case for / we arrive to the re- 
duction 

(v^A - l) [a (l - V2\) {(3 - FiOy - v^AF') - F" = 0, 

x+t(l-\/2 



where ( = — _ ^ Although this equation cannot be analytically solved 

it has for A = 0^ the special solution 



aC- 
Using it in combination with (36) we arrive at the invariant solution 



a{t + xY 

Finally, using the equivalence transformations and the boundary condition 
(4a) we obtain the similarity solution for Eq. (1) with f{() = a{( — f3)^ 

U{X,t) = P 7T 

^ a(t(a2-2r) + 21og|x|)' 

with 

H{t) = ^e^(2r-.2), 

and 

R{t) =P- ^""^ 



alog^l^l 
Case 4.2.2. /(() = -l(a + /?() (5 + log|a + /3C|) , /3,7 t^ 



■^Actually for A = its general solution can be given but in implicit form containing 
transcendental functions, hence not suitable for our analysis. 



Let the arbitrary element of the Lie algebra spanned by (25) and (26) be 
X = C1X1+C2X2+C3X3+C4X4. Using (15) and the equivalence transformations 
the first condition turns to the ODE: 

where Ci 7^ 0. The solution of (37) is 

H{t) = ^e5(2--^)*+l-(^*+^-^*), (38) 

where fx = ^^^^, A = v^(l + 2—) and A the constant of integration. Using 
the above solution, conditions (16) and (24) we get the ODE 

e*^ (k + 27/i (V2a + (77 At + 27 log A)) {a + f3R) 

+ 2e2^*7V^^(" + PR) - 8/3aR' = 0. (39) 
Its solution is 

R^t) = -^ + Be- ^ ^^i^-^ -^'" (40) 

where k = 8-^ and B the constant of integration. 

Having found the functions H, R admitted by the symmetries we proceed 
with the reduction of the Eq. (22). From the invariant surface condition (12) 
the invariant solution is 

rft + x-^t- -^e'A . (41) 

Substituting (41) to (22) for this particular case for / we arrive at the re- 
duction 

F(C) (27(25 + C) - 1 + v^A + 47(log/3 + logF(C))) + 2 (v^A - 2) F' 

- AF" = 0, 

where ( = t+x — -j=t — ^e*'''. Although for this equation the general solutions 
cannot be found in a closed form, two particular solutions can be obtained, 
one when A 7^ and one when A = 0. Each of them is considered separately: 
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F(C) = Aie-5f-i°s/^ (42) 

with Ai = e~^ . Using it in combination with (41) we arrive at the 
invariant solution 

u{x, t) = —— I Aie »T —a 

Again, using the equivalence transformations and the boundary condi- 
tion (4a) we arrive to the similarity solution for Eq. (1) with f{Q = 

-l(a + PC){5 + \og\a + m 

A e^*(CT(KCT+7M(2v^-2A-7Me^*+27<Tt))-4r72pt) 

u{x,t) = ——e ^^ Ixla^*^ — — 

with 

and 

^ eT*(<T(7M(2V^-2A + 27At + 7MeT*)+KCT)+472fjlog^) 

i?(t) = ^e ^^ - -. 

A = 

F(C) = A2e-^^^'+^^ (43) 

(7(2-4i) + (l+2C)2-47log/3) 

where C is a constant and A2 = e ^ . Using it in 

combination with (41) we arrive to the invariant solution 

M(a:,t) = e ^/' I Aae^V ^ I ' 

Again, using the equivalence transformations and the boundary con- 
dition (4a) we obtain the similarity solution for Eq. (1) with f{() = 

-l(a + ^C)iS + \og\a + PC\) 
u{x,t) = 

cr^{KtT-4s/2iiiC)e*''+l(r7'^-2A{a(',at + 2s/2(l + 2C))-2rtt)t+4j'^ log2|a;| 

A2e ^^ X 
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with 

H{t) = ^e5(2--^)*+|-Me-* 

and 



5. Conclusion 

In the present paper a generalization of the celebrated Black-Scholes- 
Merton equation (2) was proposed and studied under the prism of the mod- 
ern group analysis or symmetry method. To that end, we harnessed the 
advantage that the equivalence transformations offer when studying classes 
of differential equations, the knowledge of the best representative for this 
class of equations. This fact substantially simplifies the task of classifying it 
and obtaining its point symmetries. 

Through this classification interesting cases, from the point of symme- 
tries, arise. Nonlinear equations in general have few or no symmetries so 
cases that augment the set of symmetries at disposal are like an oasis in the 
desert. Quite commonly a dynamical system possessing an ample number 
of symmetries is more probable to relate with a physical system or model 
a more realistic process. Furthermore, in the case that we wish to study a 
boundary problem, because of the fact that not all of the symmetries admit 
the boundary and its condition, some of the symmetries will be excluded. 
Hence the bigger the set of symmetries the bigger the probability that some 
will survive the scrutiny of the boundary conditions and give an invariant 
solution for the problem in its entirety. 

For the equation studied here both sides of this fact were revealed. The 
terminal condition was too strict and gave only trivial solutions. On the 
other hand the boundary condition imposed for the barrier option allowed us 
to obtain non trivial invariant solutions, undoubtedly, the arbitrary functions 
involved in the boundary condition helped in that direction. 

The insight provided through the above symmetry analysis might prove 
practical to anyone looking for a more realistic economic model without de- 
parting from the reasoning behind the Black-Scholes-Merton equation that 
made it such a successful model on the first place. Moreover, when one 
studies more exotic kinds of options that gain ground in the Asian markets 
that in turn play an ever increasing role in the world market. Last but not 
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least although the Black-Scholes-Merton model is a standard way to price 
traditional options it encounters difficulties with exotic ones. The nonlinear 
variants of the traditional model given here, along with the found analytical 
solutions, might turn the table in that respect. We leave to the interested 
reader the possible economical interpretation and use of the obtained results. 
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